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Abstract 

m : 

Classical considerations suggest that the probability distribution P(R), where R is 
the ratio of neutral pions to total pions emitted from a disoriented chiral condensate 
C^-i, (which has been hypothesized to form in heavy ion reactions) is R 1 / 2 /2. Quantum 

mechanical isospin correlations between the condensate and the remainder of the sys- 
tern can alter this. Moments of the P(R) distribution can be expressed in terms of 

•rH , 

^ ■ expectation values of (I 2 /N 2 ) m where / is the isospin carried by the condensate, N is 

the number of pions emitted and m is an arbitrary integer. We find that the probability 
distribution is very similar to the classical distribution for 0.1 < R < 0.9 unless the 
isospin carried by the condensate is very large. 
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Recently, there has been considerable attention given to the possible formation of a 
disoriented chiral condensate (D%C) in high energy collisions of heavy ions |I|, 0, 0, f|, || 
H 0' & HI- The basic idea is rather simple. During the collision a thermalized region at 
high temperature may form. If the temperature is above the chiral restoration temperature, 
the chiral order parameters (qq) and (iq^^rq) will be zero in this region. If the system 
subsequently cools rapidly, the chiral-restored state becomes unstable. The notion of the 
disoriented chiral condensate concerns the response of the system to this instability. Clearly 
the system will eventually relax to its vacuum value in which (qq) = (qq) va , c ~ —m%f%/m q , 
and (iq r )z>Tq)= with the energy being radiated away by hadron emission. 

One scenario for this relaxation involves two time scales: a relatively quick scale in which 
a large region settles into a disoriented chiral condensate and then a slower scale in which this 
D%C relaxes to the physical vacuum. The D%C is a region which corresponds to a chirally 
rotated vacuum — i.e., a region in which the magnitude of the chiral symmetry breaking 
is the same as in the vacuum but in which the symmetry breaking is not aligned in the 
(qq) direction. The plausibility of this scenario stems from the following fact: apart from 
explicit chiral symmetry breaking (which is proportional to the current quark masses and, 
thus, small) and surface terms associated with the mismatch of the chiral direction in the 
condensate and the true vacuum (which for a large spatial volume of condensate constitutes 
a small fraction of the energy driving the system) all directions in chiral space are equally 
likely. The idea is that these effects which give a preference for "falling" in the (qq) direction 
are small perturbations on the scale of the dynamical instability and hence the system is 
almost as likely to fall in an arbitrary direction in chiral space. Of course, it is clear that 
if a region of D%C forms it will eventually decay to the physical vacuum due to the quark 
mass and surface effects. 

Whether or not this scenario is, in fact, viable is a subject of current interest |8], [5], ffi, U \. 
The key question is whether or not large regions of D%C form. Unfortunately, the detailed 
dynamics of heavy ion physics is far beyond what can be calculated directly from QCD and 
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thus, studies of this kind are, of necessity based on simplified models such as the linear sigma 



model 112 



If a large volume of D%C is formed, it will presumably decay via the emission of hadrons 
with a momentum of order V' 1 ^ 3 where V is the volume of the region with a D%C. Being the 
lightest and able to carry off chirality pions are the most likely candidates for these hadrons. 
Thus, one expects an anomalously large number of low momentum pions in the rest frame 
of D%C. A second signature concerns the relative distribution of neutral pions . 

If the low momentum pions were emitted statistically, with the pions uncorrelated with 
each other, then one would expect that the ratio R = Nq/N, of low momentum neutral pions 
to the total number of low momentum pions should have a probability distribution given by 
a binomial distribution. As N — > oo the distribution becomes sharply peaked at the average 
value, (R) = 1/3 with a variance (R 2 ) — (R) 2 = ^ which goes to zero. 

It is often asserted [0, ffl, B|, on the basis of a simple classical argument, that if the low 
mass pions were to come entirely from a D^C, then the following nonstatistical probability 
distribution emerges: 

The two are so radically different that it should be easy to tell which describes the pion 
momentum distribution better. 

The derivation of eq. ([TJ) is essentially classical |J |2|. The question we wish to address 
here is whether or not this simple form for P(R) survives a quantum treatment. It has been 
pointed out previously || [7| that if the low momentum pions in the central rapidity region 
produced during the collision from the D%C are in a quantum state of zero isospin, then the 
distribution in eq. (|l|) will be automatically reproduced in the large N limit. However, there 
is no guarantee that the low momentum pions will carry zero isospin. Indeed, in a typical 
ultrarelativistic heavy ion reaction one would expect that even if low momentum pions are 
created in large number the vast majority of particles are "high momentum" pions and other 
energetic particles. The total isospin of the system will be carried by both low momentum 



pions and the high momentum particles. While in the absence of photon emission the sum 
of these two isospins must equal the total isospin of the system, a priori there is no reason 
to assume that the low momentum pions must be in zero or any other pure isospin state. In 
order to determine how much isospin will be carried by the low momentum pions one requires 
a reliable dynamical theory of the collision process leading to emission of low momentum 
pions. Such a calculation is well beyond the current state of the art. Instead, we will simply 
cast the final state in a useful form and draw phenomenological conclusions. 

Before proceeding with our calculation it is probably useful to discuss briefly the general 
question of how the notion of a disoriented chiral condensate can be made sensible in a 
quantum mechanical context. One obvious point is that the condensed state is not the entire 
quantum mechanical state. Apart from the condensate, if formed, one expects a large number 
of statistically emitted "high energy" particles in a relativistic heavy ion reaction. Thus, the 
simplest possible description of the state would be as a product of a condensed state and a 
state describing the rest of the system. In fact, however, this is not viable. In general, if 
the chiral condensate is disoriented (i.e., not pointing in the a direction in chiral space) it 
points in some direction in isospin space, breaking isospin symmetry. We know, however, 
that the true quantum state does not break isospin (if we neglect quark mass differences 
and electromagnetic effects) . Thus, the way to incorporate the classical notion of the 
disoriented chiral condensate into a quantum picture is as follows: the physical state, to good 
approximation, may be regarded as a superposition of product states (between a condensed 
state pointing in one chiral direction and the rest of the system). This superposition involves 
a sum over states with the condensate pointing in every chiral direction with some weighting 
function. This notion of obtaining a quantum mechanically meaningful state from a mean 
field state which breaks a symmetry is quite standard. For example, in nuclear structure 
physics one may encounter a Hartree-Fock solution which is deformed, i.e., which break 
rotational symmetry. The physical states, with good quantum numbers, are obtained by 
projection which involves superposing deformed states pointing in different directions. 
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At this point, it is convenient to make a few simplifications. First, we consider collisions of 
two heavy ions of zero isospin. Second, we assume that the total isospin is conserved during 
the collision, i.e., we ignore radiative processes. Third, we assume that in such collisions 
pions are produced with essentially zero momentum in the rest frame of the condensate. 
This is possible only if the region occupied by these pions is very large. We consider events 
in which the number of zero momentum pions is large assuming that such events occur during 
heavy ion collisions. 

The concept of a disoriented chiral condensate aligned in a nontrivial isospin direction is 
meaningful only if essentially all of the pions in the condensate will be aligned in the same 
isospin direction in the condensate. This trivially implies that the pions in the condensate 
are in a completely symmetric state in isospin and by Bose symmetry they must be in a 
symmetric spatial state. All of the physical consequences discussed in this letter depend 
essentially on the low momentum pions being in an isospin symmetric state. We should note 
that the signatures discussed in this letter can tell us only whether a state of low momentum 
pions is symmetric in isospin. We have no way to determine whether such a state has its 
origin in a disoriented chiral condensate or via some other (as yet unknown) mechanism. We 
will refer to any state containing only low momentum pions in an isospin symmetric state 
as a condensate regardless of the actual mechanism by which it is formed. 

We will study moments of the probability distribution because it is possible to obtain 
compact analytical expressions for these quantities in the limit of large N. Calculation of 
P(R) as a function of R involves differences of large numbers and becomes rather unwieldy 
for N > 100. We will show only a few illustrative results for N = 100. 

To begin, we group the particles in the final state into two classes, low momentum pions 
and everything else. Assuming that the low momentum pions are all zero momentum we 
see that one completely labels the state by three quantum numbers which can be either the 
numbers of ^ and 7r°, or, the total number of pions N, the total isospin I c , and the third 
component of the isospin i 3 . We write the most general final state with zero isospin in the 
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form: 

\^> I = Q )= zl C x,ic f7rr-^ r\ N i /c ' *3>iow|-0jv; h, -«3)hi g h, (2) 

N,I c ,i 3 V U c + 1 

where Cn,i c is a numerical coefficient, the subscript "low" labels the states of the condensate 
containing low momentum pions. The notation \ipN] Ic, — * 3) high denotes the state containing 
no low momentum pions and carrying isospin I c and third component of isospin — i^; the 
subscript N merely indicates that it is the state coupling to a state of N low momentum 
pions. If all states are normalized one has the constraint that 

£icv/ = i- (3) 

N,I C 

Quantum mechanical correlations between the isospin carried by the low momentum pions 
and the rest of the system are exhibited by eq. (H). 

If all the low momentum pions are zero momentum, the low momentum pion states 
\N, I c ,i 3 )\ ow may be generated from vacuum by application of the creation operators for 
the three types of pions. Here, it is convenient to work in a cartesian basis for the isospin 
and thus the creation operators are a\ with i = 1,2,3. Given the assumption of only zero 
momentum pions, the normalized low momentum states of fixed N, I c and 23 are unique 

\N, J c ,2 3 )iow = / dft Y /c , i3 (ft) (N\)-^ 2 (n(n) ■ aY\vac) (4) 



where n(Q) is a unit vector in the Q direction and Y/ ci j 3 (f2) is the spherical harmonic function. 

Consider an operator, 0\ ow which acts only in the space of the low momentum pions, 
i.e., an operator whose matrix elements do not depend on \ipN] I c , — ^high at all. Since, the 
operator is unity in the "high" space, the expectation value of such an operator is simple: 

(Ciow) = £ |CW,/ c | 2 C'iow j (5) 

N,I C 

where 

O low (I c ) = £^y^(Ar,/ c ,z 3 |aow|iV,/ c ,23>iow (6) 

h c 

= E oJTl ^ /c ' *3l<?iow SCalar ^liV, I c , *3>low • 

h c 
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Table 1: The coefficients b™ defined in eq. © . 

Note that the number of low momentum pions N, the number of low momentum neutral 

pions iVo, and the ratio R = N /N are all examples of operators of the form of 0\ ow . 

For 0\ ow = R m , the problem of calculating the moments of the distribution is reduced to 

evaluating matrix elements R m . The general result is given below: 

~ [ ^ 2] xM f-lV+'m! // C (/ C + 1)Y 

lim R m = V V < 

N ^°° T^q J^o\ 2 2 H\j\(m - 2i)\(i - j)\[2j + 2(m - i) + 1} \ N 2 J 

It is straightforward to obtain expressions for (R m ) in the limit iV —>■ oo by using eqs 
and (0) and the obvious definition, (I 2 /N 2 ) = J2n,i c \Cn,i c \ 2 I c {I c + l)/N 2 . Writing 



(7) 



[to/2] 



lim (RT) = E b i 



(8) 



N->oo ' ' ' \ AT2i / 

the results for the coefficients b™ for m = 1 through 6 are given in Table 0. We note that 
these results agree with the classical result when (j^) 1. 

One sees from Table [I] that when (jjs) increases from to 0.1 that moments of order 
4 and higher change by as much as ~ 10%. Thus an experimental determination of high 
moments with sufficiently good accuracy can provide a measure of the isospin content of the 
condensate. 

Fortunately, there is a much more promising experimental way of relating (j^) to ob- 
servable quantities. It has been noted by Greiner, Gong and Miiller |7j that, if the low 
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momentum pions come from an I c = condensate then in every event the number of low 
momentum 7r + 's must exactly equal the number of 7r~'s. Conversely, difference between the 
two numbers provides a direct measure of the isospin content: 

, /(N + -N) 2i \ 1 / i 2i \ 

& \ — — / = 2?n \ jvs/ • (9) 

It is easy to derive the relation by noting that iV + — AL = i% and using eqs. @ and (|6|). 
Thus a measurement of ( ( - Af+ ~^"- > ) is the most direct and natural way to measure the isospin 
content of the condensate. In particular, the variance is given by 

The result in eq. ( |10|) is in a real sense model independent. It does, however, depend 
rather strongly on two underlying assumptions. The first is the assumption that a single 
large region of D%C is formed. If multiple regions form with different orientations the results 
will certainly change, and the variance in the R distribution will certainly be reduced. One 
might hope to use eq. (JTOj ) to experimentally test whether a single region of D%C is, in fact 
produced. The second important assumption is that the total system (or at least the part 
of the system which ultimately contributes) is isospin zero. This will certainly be violated 
to a certain degree. The initial state may well not be isospin zero. More importantly, 
photons will be emitted during the heavy ion reaction which violate isospin conservation. 
It is essential to estimate how large an effect these isospin violating effects will have on the 
final results. Unfortunately, estimates of the corrections to eq. (]T0|) due to isospin violation 
will, of necessity be model dependent. 

We have discussed the effect of quantum mechanical isospin correlations on moments of 
the P(R) distribution. It should be remembered, of course, that the moments are derived 
quantities — the primary quanity is P(R) itself. We would like to understand the effect of 
isospin correlations on P(R). We write the probability distribution as 

P(R) = £ \C NJ fP Ic (R) , (11) 



where the probability of finding N neutral pions in the state \N, I c , ^low is given by 

PiXR) = — ^— EK^ J c ^3|iv,/ c ,iVo)iow| 2 • (12) 



Following Horn and Silver [13], this probability is evaluated to be 



1 v y ( N-N Q + h \ y ( N-N Q -iz \ y {2I c +mN-I c )\\ 
M -2Ul^ #o! l 2 r\ 2 J' (N + I c + 1)\\ X 

/ \ 2 

2^c — is-No I ( — l) fc 2 -2 ' 



^ fa^CT I? W(/c + - , 3 - 2*)! ( ^V- 2fc )! ' (13) 

where the sums are over all integers is, j and k such that the factorials can be defined. 

We see clearly that, the quantum mechanical isospin correlations will also affect the 
probability distribution. Needless to say, we cannot calculate the distribution without a 
knowledge of the coefficients Cn,i c , which, in turn, requires a knowledge of the reaction 
mechanism. However, we can learn something by examining Pj c (R) for a range of values of 
I c . These probability distributions have sharp odd-even variations. Such variations will never 
be seen in an experiment where one measures P(R)- The sum over I c in eq. (|TT|) smoothes 
out the variations. In addition, there is bound to be additional averaging over neighboring No 
and N because of the difficulty of counting the pions exactly. We suppress these oscillations 
partially by averaging over two neighboring numbers of neutral pions. We believe that the 
procedure improves markedly the appearance of the graphs without affecting the conclusions 
we draw. In Fig.l we show these partially averaged Pj c (R) vs. R for I c = 0, 10, 20, 30 and 
60 for N = 100. The distributions for J c = 10 and 20 are not distinguishable from the 



I c = distribution which, one may recall, agrees with the classical distribution [14]. The 
distribution for I c = 30 is barely distinguishable. These graphs tell us that the classical result 
of eq.(ffl) is robust (except at the endpoint regions near R = and R — 1) and will describe 
the actual distribution P(R) unless there is massive isospin mixing. Conversely, agreement 
with the classical result will tell us little about the isospin content of the condensate. It will 
establish that the pions in the condensate are in fully isosymmetric state and, hence, in a 
fully space symmetric state. Information about the isospin admixtures is best obtained by 
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studying the charge asymmetry distribution. 
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Figure Captions 



Figure 1. P Ic (R) vs. R for I c = 0, 10, 20, 30 and 60 for N=100. The graphs for I c = 10 
and 20 are not distinguishable from that for I c = and are not labeled. Other graphs are 
labeled by their isospins. 
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